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ABSTRACT 

Background Recent studies have reported associations 
between particulate air pollution and daily mortality rates. 

Population-based, cross-sectional studies of metropolitan 
areas in the United States have also found associations 
between particulate air pollution and annual mortality rates, 
but these studies have been criticized, in part because they 
did not directly control for cigarette smoking and other 
health risks. 

Methods In this prospective cohort study, we estimated the 
effects of air pollution on mortality, while controlling for 
individual risk factors. Survival analysis, including Cox 
proportional-hazards regression modeling, was conducted 
with data from a 14-to-16-year mortality follow-up of 8111 
adults in six U.S. cities. 

Results Mortality rates were most strongly associated with cigarette smoking. After 
adjusting for smoking and other risk factors, we observed statistically significant and robust 
associations between air pollution and mortality. The adjusted mortality-rate ratio for the 
most polluted of the cities as compared with the least polluted was 1.26 (95 percent 
confidence interval, 1.08 to 1.47). Air pollution was positively associated with death from 
lung cancer and cardiopulmonary disease but not with death from other causes considered 
together. Mortality was most strongly associated with air pollution with fine particulates, 
including sulfates. <■-- 
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Conclusions Although the effects of other, unmeasured risk factors cannot be excluded with 
certainty, these results suggest that fine-particulate air pollution, or a more complex 
pollution mixture associated with fine particulate matter, contributes to excess mortality in 
certain U.S. cities. 


Several cross-sectional investigations have found associations between mortality rates and particulate 
air pollution in U.S. metropolitan areas 1,2 ’-. A recent study reported associations between infant 
mortality and particulate air pollution in the Czech Republic 4 . These studies have often been 
criticized because they did not control directly for cigarette smoking or other covariates. Recent daily 
time-series studies, which are likely to be free of confounding by individual characteristics, have 
reported associations between daily mortality rates and changes in air pollution, specifically 
particulate pollution, in London 5 and in several cities in the United States^- , -’-- , ~’--. 


Particulate air pollution is a mixture of solid particles and liquid droplets that vary in size, 
composition, and origin. Because only very small particles can be inhaled into the lungs, U.S. national 
health standards for the quality of ambient air are based on the mass concentration of "inhalable 
particles,” defined to include particles with an aerodynamic diameter of less than 10 microm- \ Fine- 
particulate air pollution includes particles with an aerodynamic diameter equal to or below 2.5 
microm. Whereas larger particles are derived chiefly from soil and other crustal materials, fine 
particles are derived primarily from the combustion of fossil fuels in transportation, manufacturing, 
and power generation. Fine-particulate pollution typically contains a mixture of particles including 
soot, acid condensates, and sulfate and nitrate particles. Fine particles are thought to pose a 
particularly great risk to health because they are more likely to be toxic than larger particles and can 
be breathed more deeply into the lungs 1 - . 

In this study, a well-characterized cohort of adults participating in the Harvard Six Cities Study of the 
health effects of air pollution was followed prospectively, beginning in 1974—. The objective of this 
study was to estimate the effects of air pollution on mortality, with control for individual smoking 
status, sex, age, and other risk factors. 

Methods 

Study Population 

We selected random samples of adults from six communities—: Watertown, Massachusetts (where 
study enrollment was conducted in 1974); Harriman, Tennessee, including Kingston (1975); specific 
census tracts of St. Louis (1975); Steubenville, Ohio (1976); Portage, Wisconsin, including Wyocena 
and Pardeeville (1976); and Topeka, Kansas (1977). The sample was restricted to the 8111 white 
subjects who were 25 through 74 years of age at enrollment, had undergone spirometric testing, and 
had completed a standardized questionnaire. The questionnaire included questions about age, sex, 
weight, height, education level, complete smoking history, occupational exposures, and medical 
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history. 

Informational letters and postage-paid return postcards including a question on vital status were 
mailed to the subjects annually. The vital status of the subjects who did not respond was determined 
by questioning family members, friends, or neighbors. In addition, we searched the National Death 
Index -- for the years 1979 through 1989. Death certificates were obtained for 1401 of the 1430 
subjects who had died (98 percent); the causes of death were coded according to the International 
Classification of Diseases, 9th Revision (ICD-9) by an independent certified nosologist who was 
blinded both to pollution levels and to the study design and objectives. The ending date of the study 
for each city was March or June of 1991, depending on the date of the last follow-up contact; the total 
duration of follow-up was 14 to 16 years (111,076 person-years). 

For subjects who died, survival times were calculated by subtracting the date of enrollment from the 
exact date of death. For surviving participants who were not lost to follow-up, censored survival 
times were defined as the date of the end of the study minus the enrollment date. For those who were 
lost to follow-up before the period covered by our National Death Index search (i.e., before 1979), 
censored survival times were estimated by subtracting the enrollment date either from the date of the 
last follow-up contact plus six months or from the first day of the National Death Index search period 
(January 1,1979), whichever came first. For those who were lost to follow-up after the National 
Death Index search period (i.e., after 1989), censored survival times were estimated by subtracting the 
enrollment date either from the date of the last follow-up contact plus six months or from the last day 
of the study period, whichever came first. For those who were lost to follow-up during the period 
covered by the National Death Index search, the censored survival times were estimated by 
subtracting the date of enrollment from the last date in the search period (December 31, 1989). 

Air-Pollution Data 

As part of the original study design, ambient (outdoor) concentrations of total suspended particulate 
matter, sulfur dioxide, ozone, and suspended sulfates were measured in each community at a centrally 
located air-monitoring station--. Size-selective aerosol samplers were placed at these sites in the late 
1970s; data were collected for two classes of particle: fine particles (aerodynamic diameter <2.5 
microm) and inhalable particles (aerodynamic diameter, <15 microm before 1984 and <10 microm 
starting in 1984). In the mid-1980s, supplemental 24-hour integrated sampling of aerosol acidity by 
the measurement of hydrogen ion concentrations 1 - was conducted for approximately one year in each 
city. Mean pollution levels for each pollutant were calculated for periods that were consistent and 
comparable among the six cities. 

Statistical Analysis 

Life-table survival probabilities for each year of follow-up were estimated for each city, and 

1 Q 

differences between city-specific mortality rates were assessed with a log-rank test- . We estimated 
adjusted mortality-rate ratios for air pollution by simultaneously adjusting for other risk factors in 
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Cox proportion a I-hazards regression models In these models the subjects were stratified 

according to sex and five-year age groups, and each sex-age group had its own base-line hazard. Each 
model also included indicator variables for current or former smokers, the number of pack-years of 
smoking (evaluated separately for current and former smokers), an indicator variable for less than a 
high-school education, and body-mass index (defined as the weight in kilograms divided by the 
square of the height in meters). 

Two approaches were used to evaluate the effects of air pollution in the Cox proportional-hazards 
models. First, indicator variables for the city of residence were included, with Portage, Wisconsin, the 
city with the lowest levels of particulate air pollution, as the reference category. Adjusted mortality- 
rate ratios for each of the six cities were then compared graphically with the mean pollution levels in 
those cities. Next, adjusted mortality-rate ratios were estimated by including city-specific pollution 
levels directly in the Cox proportional-hazards models. Adjusted rate ratios were calculated and 
reported for a difference in air pollution equal to that between the city with the highest levels of air 
pollution and the city with the lowest levels -- that is, the adjusted rate ratios across the range of 
exposure for each pollutant among the six cities. 

Analyses were conducted to evaluate the robustness of the models and the possibility of residual 
confounding. Models were estimated after the data were separated according to the subjects' smoking 
status, sex, and occupational exposure to dust, gases, or fumes. The effect of the inclusion of different 
covariates on the estimated effect of pollution was evaluated. Models were also estimated after the 
exclusion of subjects who had been treated for high blood pressure within 10 years of enrollment in 
the study and subjects who had ever been told by a doctor that they had diabetes, had glucose in their 
urine, or had too much glucose in their blood. We also used a variety of approaches to estimate 
censored survival times. 

Mortality-rate ratios from the Cox proportional-hazards models (with adjustment for cigarette 
smoking, education, and body-mass index) were estimated separately for the following cause-of-death 
categories: cardiopulmonary (ICD-9 codes 400 through 440 and 485 through 496), lung cancer (162), 
and all others. For each cause-of-death category, data on subjects whose deaths were not in that 
specific category were censored at the time of death. 

Results 

Characteristics of the Cohort and Air-Pollution Data 

The characteristics of the cohort and the values for air-pollution measures are summarized in Ta ble 1 . 

For all measures of air pollution except the ozone level and aerosol acidity, ambient concentrations 
were highest in Steubenville and lowest in Portage or Topeka. The mean acidity of the aerosol was 
highest in Harriman, but second-highest in Steubenville. The mean ozone concentrations were highest 
in Portage and Topeka. The concentrations of total particles declined during the study period, 
especially in Steubenville and St. Louis; the annual average concentrations of fine and sulfate 
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particles varied relatively little during the study period (Figure 1). Crude mortality rates (Table 1) and 
survival curves (Figure 2) both show that mortality was highest in Steubenville and St. Louis and 
lowest in Portage and Topeka. Differences in the probability of survival among the cities were 
statistically significant (PcO.OOl). 


( View this table: Table 1. Characteristics of l he Study Population and Mean Air-Poll ill ion 
j [in t his w indo w ] Levels in Six Cities, 
jfi n a ne w windo w] 

1 



Figure 1. Annual Average Concentrations of Total Particles, Fine 
Particles, and Sulfate Particles in the Six Cilies. 



Figure 2, Crude Probability of Survival in the Six Cities; 
According to Years of Follow -up. 
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Adjusted Mortality Rates 

On the basis of the proportional-hazards model, mortality was most strongly associated with cigarette 
smoking ( Table 2). Increased mortality was also associated with having less than a high-school 
education and with increased body-mass index (the latter was especially true for women). After 
simultaneous adjustment for these other risk factors, the differences in mortality among the six cities 
remained significant. 


j; View this table: Table 2. Adjusted Mortality-Rate Ratios Estimated from Cox Proportional - 
1 [in this window] Hazards Models. 

Lin a new window] j A 


City-specific mortality rates, adjusted for a variety of health risk factors, were associated with the 
average levels of air pollutants in the cities ( Fi g u r e 3 ). The small differences in ozone levels among 
the cities (Table 1) limited the power of the study to detect associations between mortality and ozone 
levels. Mortality was more strongly associated with the levels of inhalable, fine, and sulfate particles 
than with the levels of total suspended particles, the sulfur dioxide levels, the nitrogen dioxide levels, 
or the acidity of the aerosol. 



(25K): 

[in this win dow ] 
fin a new window 1 


Figure 3. Estimated Adjusted Mortality-Rate Ratios and Pollution 
Levels in the Six Cities. 

Mean values are shown for the measures of air pollution. P denotes 
Portage. Wisconsin; T Topeka, Kansas; W Watertown, 
Massachusetts; L St. Louis; 11 Harriman, Tennessee; and S 
Steubenville, Ohio. 


When the mean concentrations of each pollutant were included individually in the proportional- 
hazards model, we found significant associations between mortality and inhalable, fine, or sulfate 
particles (P<0.005). For a difference in the air-pollution level equal to that between the most polluted 
city and the least polluted city and with inhalable particles (range, 18.2 to 46.5 jig per cubic meter), 
fine particles (range, 11.0 to 29.6 g.g per cubic meter), and sulfate particles (range, 4.8 to 12.8 /.rg per 
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cubic meter) used as indicators of air pollution, the adjusted rate ratios were nearly equal at 1.27 (95 
percent confidence interval, 1.08 to l .48), 1.26 (95 percent confidence interval, 1.08 to 1.4-7), and 
1.26 (95 percent confidence interval, 1.08 to 1.47), respectively. 

Sensitivity 

Estimates of the association between mortality and fine-particle pollution among subjects with 
different smoking status and among men and women (Table 3) showed only small and nonsignificant 
differences between subgroups. Associations with air pollution were somewhat stronger among 
subjects with occupational exposure to dust, gases, or fumes (T able 3). However, positive 
associations between mortality and air-pollution levels were observed in all subgroups defined by 
occupational exposure and sex, and differences among the subgroups were not statistically 
significant. 


View this table: Table 3. Adjusted Mortality-Rate, Ratios for the Most Polluted and Least 
f in this w indow 1 Polluted Cities Studied, According to Smoking Status, Sex, and 
1 i n a new window] Occupational Exposure, with Fine Particles Used as the Indicator of Air 
Pollution. 


Although cigarette smoking and other risk factors were associated with mortality, our estimates of 
pollution-related mortality were not significantly affected by the inclusion or exclusion of these 
variables in the models (Table 4). The estimated association of air pollution and mortality was 
unchanged when subjects who had been treated for high blood pressure or subjects with diabetes were 
excluded from the analysis (T able 4 ). When censored survival times were recalculated as the date of 
the last follow-up contact minus the enrollment date, or when the analysis was restricted to data on 
deaths in 1979 through 1989 (the years of the National Death Index searches), no appreciable 
differences in the estimated association between air pollution and mortality were observed. 


View this table: Table 4. Pstimaied Mortality-Rate Ratios for the Most Polluted City as 
1 fin this window 1 i Compared with the Least Polluted City, with Fine Particles Used as the. 
1 in a new' window I Indicator of Air Pollution, in Selected Models. 


Causes of Death 

The estimated effects of air pollution on mortality varied among causes of death (Ta bl e 5). For 
comparison, rate ratios were estimated for current smokers and for former smokers with 
approximately the average number of pack-years of smoking at enrollment (Table 5). Smoking was 
most strongly associated with mortality due to lung cancer, significantly associated with mortality due 
to cardiopulmonary disease, but not associated with mortality from all other causes. Similarly, air 
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pollution was positively associated with mortality due to lung cancer and cardiopulmonary disease 
but not with mortality from all other causes. Only 98 deaths were coded on the death certificates as 
due to nonmalignant respiratory disease (ICD-9 codes 485 through 496), as compared with 646 deaths 
due to cardiovascular disease (codes 400 through 440). An analysis restricted to deaths from 
nonmalignant respiratory disease produced unstable and statistically nonsignificant estimates of the 
association with air pollution. When mortality from all causes was considered, or when deaths due to 
cardiovascular and respiratory diseases were grouped together, the effects of air pollution were 
consistent and the association was robust. 


| View this table: Table 5. Adjusted Mortality-Rate Ratios for Current and Former Cigarette 
I fin this window! Smokers and for the Most Polluted City as Compared with the Least 
ijtn a new window 1 Polluted, According to Cause of Death. 


Discussion 

In this prospective cohort study, the mortality rate, adjusted for other health risk factors, was 
associated with the level of air pollution. Mortality was more strongly associated with the levels of 
fine, inhalable, and sulfate particles than with the levels of total particulate pollution, aerosol acidity, 
sulfur dioxide, or nitrogen dioxide. As with all other epidemiologic studies, it is possible that the 
observed association was due to confounding — that is, that it resulted from a risk factor that was 
correlated with both exposure and mortality. Potential confounders of the effects of air pollution 
include cigarette smoking and occupational exposure to pollutants. In our study, however, the 
association of air pollution with mortality was observed even after we directly controlled for 
individual differences in other risk factors, including age, sex, cigarette smoking, education level, 
body-mass index, and occupational exposure. 

The estimated effect of air pollution on mortality was not altered by the inclusion or exclusion of 
indicator variables for other risk factors in our models. Analyses were conducted for Subgroups 
defined according to sex, smoking status, and occupational exposure. Although the effects of 
pollution were somewhat stronger among subjects occupationally exposed to dust, gases, or fumes, 
positive associations between mortality and air pollution were observed among all the smoking- 
status, occupational-exposure, and sex groups, and the differences among these subgroups were not 
statistically significant. The estimated association of pollution and mortality remained essentially 
unchanged when subjects who had been treated for high blood pressure or who had diabetes were 
excluded from the analysis. 

In our analysis, the mortality-rate ratios have been expressed in terms of the range of exposure to air 
pollutants in the six cities. When the range of exposure was used, the estimated relative rate ratios for 
inhalable, fine, and sulfate particles were nearly equal at 1.27 (95 percent confidence interval, 1.08 to 
1.48), 1.26 (95 percent confidence interval, 1.08 to 1.47), and 1.26 (95 percent confidence interval, 
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1.08 to 1.47), respectively, Because the six cities were selected as representative of the range of 
particulate air pollution in the United States, these rate ratios roughly represent the relative risk 
associated with that range. 

In this study, exposure to air pollution was estimated by monitoring outdoor air pollution at a central 
site in each of the six cities. Long-term transport and large-scale mixing of combustion products play 
a large part in establishing the levels of sulfate and fine-particulate air pollution. Therefore, 
concentrations of sulfates and fine particles are relatively uniform within each of these 
communities*"'. Furthermore, sulfate and fine-particulate air pollution penetrates indoors, resulting in 
strong correlations between indoor and outdoor concentrations^ 4 ’--’—. Thus, measurements of the 
outdoor concentrations of sulfate and fine particles may be better indicators of individual exposure 
than the other pollutants we considered. 

The associations observed in this study between air pollution and mortality are consistent with 
associations observed in recent time-series studies, including studies from three of these six 
cities^’ 6 ’ 7 ’®’-’—Because the daily time-series studies evaluated only the effect of short-term 
changes in pollution levels, whereas our study evaluated associations with long-term exposure 
(including recurring episodes of relatively high pollution), quantitative comparisons with these 
investigations are difficult to make. Nevertheless, as was found in the time-series studies, particulate 
air pollution was associated with death due to cardiopulmonary causes. In our study, in which we 
evaluated the effects of long-term exposure, lung cancer was associated with particulate air pollution; 
such an association with lung cancer was not observed in the daily time-series studies. Little or no 
association with other causes of death was evident in our study or the time-series studies. The small 
number of reported deaths due to nonmalignant respiratory disease and the potential for 
misclassification of primary causes inherent in the use of death-certificate data limited our ability to 
evaluate cause-specific mortality in more detail. 

The pollution concentrations used in our analysis represent only exposures monitored during the 
study period. Increased mortality, however, may reflect the cumulative burden of a lifetime of 
exposure. Concentrations of total particles clearly declined during the study period (especially in 
Steubenville and St. Louis), whereas concentrations of fine particles and sulfate particles were 
relatively stable. Given the lack of data on pollution levels before the study period and in view of the 
fact that the relative ranking of the cities in terms of air-pollution levels did not change during the 
study period, it is not possible to differentiate the influences of historical exposure from those of 
recent exposure. The observed association between mortality and mean exposure to fine-particulate 
and sulfate air pollution during the study period may also partially reflect exposure to air pollution 
before the study period. 

The strength of the observed association between air pollution and mortality is confirmed by previous 
observations of associations between particulate air pollution and other health end points. Elevated 
levels of particulate air pollution have been associated with declines in lung function or with 
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increases in respiratory symptoms such as cough, shortness of breath, wheezing, and asthma 


27.2 S .29.30.31.32.33.34.35.36 


. Other studies have found associations between particulate air 


attacks' 

pollution and rates of hospitalization/ 7 ” 1 -’- 9-.40.4I chronic obstructive pulmonary disease, 4 -? 
restricted activity due to illness 4 -’- 4 . 


and 


A large and growing body of literature documents the adverse health effects associated with 
particulate air pollution. Although the effects of unmeasured risk factors cannot be controlled for, in 
this prospective cohort study we observed significant effects of air pollution on mortality even when 
we controlled for sex, age, smoking status, education level, and occupational exposure to dust, gases, 
and fumes. The compatibility of the effects of air pollution on mortality in this study with those 
observed in population-based cross-sectional studies and daily time-series studies provides further 
evidence for the conclusion that exposure to air pollution contributes to excess mortality. This study, 
therefore, provides additional impetus to the development of strategies to reduce urban air pollution. 

Supported in part by grants (ES-01108 and ES-00002) from the National Institute of Environmental Health Sciences, by 
cooperative agreements (CR-811650 and CR-818090) with the Environmental Protection Agency, and by a contract (RP- 
1001) with the Electric Power Research Institute. 
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